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ABSTRACT 

Regular chains and triangular decompositions are fundamen- 
tal and well-developed tools for describing the complex so- 
lutions of polynomial systems. This paper proposes adapta- 
tions of these tools focusing on solutions of the real analogue: 
semi-algebraic systems. 

We show that any such system can be decomposed into 
finitely many regular semi-algebraic systems. We propose 
two specifications of such a decomposition and present cor- 
responding algorithms. Under some assumptions, one type 
of decomposition can be computed in singly exponential time 
w.r.t. the number of variables. We implement our algorithms 
and the experimental results illustrate their effectiveness. 

Categories and Subject Descriptors 

1.1.2 [Symbolic and Algebraic Manipulation]: Algo- 
rithms — Algebraic algorithms, Analysis of algorithms 

General Terms 

Algorithms, Experimentation, Theory 

Keywords 

regular semi-algebraic system, regular chain, triangular de- 
composition, border polynomial, fingerprint polynomial set 

1. INTRODUCTION 

Regular chains, the output of triangular decompositions 
of systems of polynomial equations, enjoy remarkable prop- 
erties. Size estimates play in their favor Q/J] and permit the 
design of modular [J3] and fast [17] methods for computing 
triangular decompositions. These features stimulate the de- 
velopment of algorithms and software for solving polynomial 
systems via triangular decompositions. 

For the fundamental case of semi-algebraic systems with 
rational number coefficients, to which this paper is devoted, 
we observe that several algorithms for studying the real so- 
lutions of such systems take advantage of the structure of a 
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regular chain. Some are specialized to isolating the real solu- 
tions of systems with finitely many complex solutions [231 1101 
[3]. Other algorithms deal with parametric polynomial sys- 
tems via real root classification (RRC) [2S] or with arbitrary 
systems via cylindrical algebraic decompositions (CAD) [9]. 

In this paper, we introduce the notion of a regular semi- 
algebraic system, which in broad terms is the "real" counter- 
part of the notion of a regular chain. Then we define two 
notions of a decomposition of a semi- algebraic system: one 
that we call lazy triangular decomposition, where the ana- 
lysis of components of strictly smaller dimension is deferred, 
and one that we call full triangular decomposition where all 
cases are worked out. These decompositions are obtained by 
combining triangular decompositions of algebraic sets over 
the complex field with a special Quantifier Elimination (QE) 
method based on RRC techniques. 

Regular semi-algebraic system. Let T be a regular chain of 
Q[a;r, . • • , x n ] for some ordering of the variables x = xi, . . . , x n . 
Let u = ui, . . . ,Ud and y = yi, . . . , y n -d designate respec- 
tively the variables of x that are free and algebraic w.r.t. T. 
Let P C Q[x] be finite such that each polynomial in P is 
regular w.r.t. the saturated ideal of T. Define P > := {p > 
| p £ P}. Let Q be a quantifier-free formula of Q[x] in- 
volving only the variables of u. We say that R := [Q, T, P>] 
is a regular semi-algebraic system if: 
(i) Q defines a non-empty open semi-algebraic set S in R d , 
(ii) the regular system [T, P] specializes well at every point 

u of S (see Section [2] for this notion), 
(iii) at each point u of S, the specialized system [T(u), P(ii)>] 

has at least one real zero. 
The zero set of R, denoted by Zr(R), is defined as the set 
of points (u, y) £ R d x R n ~ d such that Q(u) is true and 
t{u, y) = 0, p(u, y) > 0, for all t € T and all p e P. 

Triangular decomposition of a semi-algebraic system. In 

Section we show that the zero set of any semi-algebraic 
system 6 can be decomposed as a finite union (possibly 
empty) of zero sets of regular semi-algebraic systems. We 
call such a decomposition a full triangular decomposition (or 
simply triangular decomposition when clear from context) of 
&, and denote by RealTriangularize an algorithm to compute 
it. The proof of our statement relies on triangular decompo- 
sitions in the sense of Lazard (see Section [2] for this notion) 
for which it is not known whether or not they can be com- 
puted in singly exponential time w.r.t. the number of vari- 
ables. Meanwhile, we are hoping to obtain an algorithm for 
decomposing semi-algebraic systems (certainly under some 
genericity assumptions) that would fit in that complexity 
class. Moreover, we observe that, in practice, full triangular 



decompositions are not always necessary and that providing 
information about the components of maximum dimension 
is often sufficient. These theoretical and practical motiva- 
tions lead us to a weaker notion of a decomposition of a 
semi-algebraic system. 

Lazy triangular decomposition of a semi-algebraic sys- 
tem. Let & — [F, JV>, P>, Hjt] (see Section for this nota- 
tion) be a semi-algebraic system of Q[x] and Z^(&) C R n be 
its zero set. Denote by d the dimension of the constructible 
set (i 6 C" f(x) = 0,g(x) # 0, for all / G F,g £ P U H}. 
A finite set of regular semi-algebraic systems E4, i = 1 • • • t 
is called a lazy triangular decomposition of 6 if 

• Ul =1 Zu(Ri) C Zk(&) holds, and 

• there exists G C Q[x] such that the real-zero set Zr(G) 
C I" contains Z R (6) \ (u' =1 Z K (_Ri)) and the complex- 
zero set V(G) C C™ either is empty or has dimension 
less than d. 

We denote by LazyRealTriangularize an algorithm comput- 
ing such a decomposition. In the implementation presented 
hereafter, LazyRealTriangularize outputs additional informa- 
tion in order to continue the computations and obtain a full 
triangular decomposition, if needed. This additional infor- 
mation appears in the form of unevaluated function calls, 
explaining the usage of the adjective lazy in this type of de- 
compositions. 

Complexity results for lazy triangular decomposition. In 

Section |4j we provide a running time estimate for computing 
a lazy triangular decomposition of the semi-algebraic sys- 
tem 6 when & has no inequations nor inequalities, (that is, 
when N> — P> = — holds) and when F generates 
a strongly equidimensional ideal of dimension d. We show 
that one can compute such a decomposition in time singly 
exponential w.r.t. n. Our estimates are not sharp and are 
just meant to reach a singly exponential bound. We rely 
on the work of J. Renagar [20] for QE. In Sections [5] and [6] 
we turn our attention to algorithms that are more suitable 
for implementation even though they rely on sub-algorithms 
with a doubly exponential running time w.r.t. d. 

A special case of quantifier elimination. By means of tri- 
angular decomposition of algebraic sets over C, triangular 
decomposition of semi-algebraic systems (both full and lazy) 
reduces to a special case of QE. In Section [SJ we implement 
this latter step via the concept of a fingerprint polynomial 
set, which is inspired by that of a discrimination polynomial 
set used for RRC in [231 |2"4] . 

Implementation and experimental results. In Section [S] 
we describe the algorithms that we have implemented for 
computing triangular decompositions (both full and lazy) of 
semi-algebraic systems. Our Maple code is written on top 
of the RegularChains library. We provide experimental data 
for two groups of well-known problems. In the first group, 
each input semi-algebraic system consists of equations only 
while the second group is a collection of QE problems. To 
illustrate the difficulty of our test problems, and only for this 
purpose, we provide timings obtained with other well-known 
polynomial system solvers which are based on algorithms 
whose running time estimates are comparable to ours. For 
this first group we use the Maple command Groebner : - 
Basis for computing lexicographical Grobner bases. For the 
second group we use a general purpose QE software: QEPCAD 
b (in its non-interactive mode) [5] . Our experimental results 
show that our LazyRealTriangularize code can solve most of 



our test problems and that it can solve more problems than 
the package it is compared to. 

We conclude this introduction by computing a triangular 
decomposition of a particular semi-algebraic system taken 
from [BJ. Consider the following question: when does p(z) = 
z 3 + az + b have a non-real root x + iy satisfying xy < 1 ? 
This problem can be expressed as (3x)(3y)[f = g — A y 7^ 
f\xy — 1 < 0], where / = Ke(p(x + iyj) = x s — 3xy 2 + ax + b 
and g = lm(p(x + iy))/y = 3a; 2 — y 2 + a. 

We call our LazyRealTriangularize command on the semi- 
algebraic system / = 0, g — 0, y ^ 0, xy — 1 < with the 
variable order y > x > b > a. Its first step is to call the 
Triangularize command of the RegularChains library on the 
algebraic system / = g = 0. We obtain one squarefree regu- 
lar chain T = \ti , £2] , where ti = g and t 2 = 8a; 3 + 2aa; — b, 
satisfying V(f,g) — V(T). The second step of LazyRe- 
alTriangularize is to check whether the polynomials defin- 
ing inequalities and inequations are regular w.r.t. the sat- 
urated ideal of T, which is the case here. The third step 
is to compute the so called border polynomial set (see Sec- 
tion [H which is B = [hi,h 2 ] with hi = 4a 3 + 27b 2 and 
h 2 = -4a 3 6 2 - 276 4 + 16a 4 +512a 2 +4096. One can check that 
the regular system [T, {y, xy — 1}] specializes well outside of 
the hypersurface /11/12 = 0. The fourth step is to compute 
the fingerprint polynomial set which yields the quantifier- 
free formula Q = hi > telling us that [Q, T, 1 — xy > 0] is 
a regular semi-algebraic system. After performing these four 
steps, (based on Algorithm [S] Section [SJ) the function call 

LazyRealTriangularize([/, g, y ^ 0, xy — 1 < 0], [y, x, b, a]) 

in our implementation returns the following: 

' [{t! = 0, t 2 = 0, 1 - xy > 0, hi > 0}] h 1 h 2 ^0 

•/.LazyRealTriangularize ([ti = 0, t 2 = 0, / = 0, 
• hi = 0, 1-xy > 0,y ^ 0],[y,x,b,a]) hi = 

"/.LazyRealTriangularize ([ti = 0, t 2 = 0, / = 0, 

hz = 0, 1 - xy > 0,y ^ 0], [y,x,b,a]) h 2 = 

The above output shows that {[Q,T, 1 — xy > 0]} forms 
a lazy triangular decomposition of the input semi-algebraic 
system. Moreover, together with the output of the recursive 
calls, one obtains a full triangular decomposition. Note that 
the cases of the two recursive calls correspond to hi = and 
h 2 = 0. Since our LazyRealTriangularize uses the Maple 
piecewise structure for formatting its output, one simply 
needs to evaluate the recursive calls with the value com- 
mand, yielding the same result as directly calling RealTrian- 
gularize 

C [{h = 0, t 2 = 0, 1 - xy > 0, hi > 0}] /ii/i 2 ^0 

[] /u=o 
' c ([{t 3 =o,t 4 = 0, 

I h 2 = 0,l-xy>0}}) h 3 ^0 h2 = (j 
■ [ [] h 3 = 

where t 3 = xy + 1, U = 2a? x - a 2 b + 32ax - 486 + 18a:fo 2 , 
h 3 = (a 2 + 48)(a 2 + 16)(a 2 + 12). 

From this output, after some simplification, one could ob- 
tain the equivalent quantifier-free formula, 4a 3 + 27& 2 > 0, 
of the original QE problem. 



2. TRIANGULAR DECOMPOSITION OF 
ALGEBRAIC SETS 

We review in the section the basic notions related to reg- 
ular chains and triangular decompositions of algebraic sets. 
Throughout this paper, let k be a field of characteristic and 
K be its algebraic closure. Let k[x] be the polynomial ring 
over k and with ordered variables x = xi < ■ ■ ■ < x n . Let 
p,q 6 k[x] be polynomials. Assume that p ^ k. Then denote 
by mvar(p), init(p), and mdeg(p) respectively the greatest 
variable appearing in p (called the main variable of p), the 
leading coefficient of p w.r.t. mvar(p) (called the initial of p) , 
and the degree of p w.r.t. mvar(p) (called the main degree 
of p); denote by der(p) the derivative of p w.r.t. mvar(p); 
denote by discrim(p) the discriminant of p w.r.t. mvar(p). 

Triangular sets. Let T C k[x] be a triangular set, that is, a 
set of non-constant polynomials with pairwise distinct main 
variables. Denote by mvar(T) the set of main variables of 
the polynomials in T. A variable v in x is called algebraic 
w.r.t. T if v £ mvar(T), otherwise it is said free w.r.t. T. 
If no confusion is possible, we shall always denote by u = 
Mi , . . . , Ud and y = yi , . . . , y m respectively the free and the 
main variables of T. Let hr be the product of the initials 
of the polynomials in T. We denote by sat(T) the saturated 
ideal of T: if T is the empty triangular set, then sat(T) 
is defined as the trivial ideal (0), otherwise it is the ideal 
(T) : fey. The quasi- component W(T) of T is defined as 
V(T) \ V(h T ). Denote W(T) = V(sat(T)) as the Zariski 
closure of W(T). 

Iterated resultant. Let p and q be two polynomials of k[x]. 
Assume q is non-constant and let v = invar (g). We de- 
fine res(p, q, v) as follows: if v does not appear in p, then 
res(p, q, v) := p; otherwise res(p, q, v) is the resultant of p 
and q w.r.t. v. Let T be a triangular set of k[x]. We define 
res(p, T) by induction: if T is empty, then res(p,T) = p; 
otherwise let v be the greatest variable appearing in T, then 
res(p, T) = res(res(p, T v , v), T<„), where T v and T<„ denote 
respectively the polynomials of T with main variables equal 
to and less than v. 

Regular chain. A triangular set T C k[x] is a regular chain 
if: either T is empty; or (letting t be the polynomial in T 
with maximum main variable), T \ {t} is a regular chain, 
and the initial of t is regular w.r.t. sat(T \ {£}). The empty 
regular chain is denoted by 0. Let H C k[x]. The pair 
[T, H] is a regular system if each polynomial in H is regular 
modulo sat(T). A regular chain T or a regular system [T, H], 
is squarefree if for all t € T, the der(t) is regular w.r.t. sat(T). 

Triangular decomposition. Let F C k[x]. Regular chains 
Ti, . . . , T e of k[x] form a triangular decomposition of V(F) if 
either: V(F) = Uf =1 W(Ti) (Kalkbrener's sense) or V(F) = 
Ui = iW(Ti) (Lazard's sense). In this paper, we denote by 
Triangularize an algorithm, such as the one of |18| . computing 
a triangular decomposition of the former kind. 

Regularization. Let p 6 k[x]. Let T be a regular chain of 
k[x]. Denote by Regularize^, T) an operation which com- 
putes a set of regular chains {Ti, . . . ,T e } such that (1) for 
each i,l < i < e, either p £ sat(T») or p is regular w.r.t. 
sat(Ti); (2) we have W(T) = W(Ti)lJ- ■ -UW(T e ), mvar(T) = 
mvar(Ti) and sat(T) C sat(Ti) for 1 < i < e. 

Good specialization [8]. Consider a squarefree regular sys- 
tem [T, H] of k[u, y] . Recall that y and u = ui, ■ ■ ■ , Wd stand 
respectively for mvar(T) and x \ y. Let z = (21, ... , Zd) be 



a point of K d . We say that [T, H] specializes well at z if: (i) 
none of the initials of the polynomials in T vanishes mod- 
ulo the ideal (zi — Ui, . . ■ , Zd — Ud)\ {H) the image of [T, H] 
modulo {zi — iti, . . . , Zd — Ud) is a squarefree regular system. 

Border polynomial | 25 | . Let [T, H] be a squarefree regular 
system of k[u, y]. Let bp be the primitive and square free 
part of the product of all res(der(i), T) and all res(h,T) for 
h £ H and t 6 T. We call bp the border polynomial of [T, H] 
and denote by BorderPolynomial (T, H) an algorithm to com- 
pute it. We call the set of irreducible factors of bp the border- 
polynomial set of [T, H]. Denote by BorderPolynomialSet(T, H) 
an algorithm to compute it. Proposition [1] follows from the 
specialization property of subresultants and states a funda- 
mental property of border polynomials. 

Proposition 1. The system [T, H] specializes well atu € 
K d if and only if the border polynomial bp{u) 7^ 0. 

3. TRIANGULAR DECOMPOSITION OF 
SEMI-ALGEBRAIC SYSTEMS 

In this section, we prove that any semi-algebraic system 
can be decomposed into finitely many regular semi-algebraic 
systems. This latter notion was defined in the introduction. 
Semi-algebraic system. Consider four finite polynomial 
subsets F = {fi,fa, ■•• , /<,}, N = {ni,n 2 , ■ ■ ■ ,n t }, P = 
{pi,P2, ■ ■ ■ ,Pr}, and H = {hi,h 2 , - ■ ■ ,h e } oi Q[xi, ...,x„]. 
Let iV> denote the set of non-negative inequalities {rii > 
0, . . . ,n t > 0}. Let P> denote the set of positive inequal- 
ities {pi > 0,...,p r > 0}. Let denote the set of in- 
equations {hi 7^ 0, . . . , he 7^ 0}. We will denote by [F, P>] 
the basic semi-algebraic system {/1 = 0, . . . , f s — 0,pi > 
0,...,p r > 0}. We denote by 6 = [F, N> , P> , H^t] the 
semi-algebraic system (SAS) which is the conjunction of the 
following conditions: /1 = 0, . . . , f s — 0, n\ > 0, . . . , nt > 0, 
pi > 0, . . . ,p r > and hi ^ 0, . . . , hi ^ 0. 

Notations on zero sets. In this paper, we use U Z" to denote 
the zero set of a polynomial system, involving equations and 
inequations, in C n and "Zr" to denote the zero set of a semi- 
algebraic system in R n . 

Pre-regular semi-algebraic system. Let [T, P] be a square- 
free regular system of Q[u, y]. Let bp be the border polyno- 
mial of [T, P]. Let B C Q[u] be a polynomial set such that bp 
divides the product of polynomials in B. We call the triple 
[B^tjT, P>] a pre-regular semi- algebraic system of Q[x]. Its 
zero set, written as Zm(Bjt, T, P>), is the set (u,y) £ R n 
such that b(u) 7^ 0, t(u,y) = 0, p(u,y) > 0, for all b 6 B, 
t G T, p 6 P. Lemma [T] and Lemma [2] are fundamental 
properties of pre-regular semi-algebraic systems. 

Lemma 1. Let & be a semi- algebraic system of Q[x]. Then 
there exists finitely many pre-regular semi-algebraic systems 
[Bi^T^P^], i = l-.-e, s.t. Zm(&) = Uf = iZ R (B^, Ti, P %> ). 

PROOF. The semi-algebraic system 6 decomposes into 
basic semi-algebraic systems, by rewriting inequality of type 
n > as: n > V n = 0. Let [F, P>] be one of those 
basic semi-algebraic systems. If F is empty, then the triple 
[P, 0,P>], is a pre-regular semi-algebraic system. If F is 
not empty, by Proposition [T] and the specifications of Tri- 
angularize and Regularize, one can compute finitely many 
squarefree regular systems [Ti,H] such that V(F)C\Z(P 1 t) = 
Uf = i (V(Tj) n Z(Bi^)) holds and where B r is the border 
polynomial set of the regular system [Ti, H]. Hence, we have 



Z R (P,P>) = U? =1 Z R (B i# ,Ti,P>), where each [P^,P,P>] 
is a pre-regular semi-algebraic system. □ 

Lemma 2. Let [B-t, T, P>] be a pre-regular semi-algebraic 
system o/Q[u, y]. Let h be the product of polynomials in B. 
The complement of the hypersurface h — in R d consists 
of finitely many open cells of dimension d. Let C be one 
of them. Then for all a £ C , the number of real zeros of 
[T(a), P>(a)] is the same. 

Proof. From Proposition[T]and recursive use of Theorem 
1 in [11] on the delineability of a polynomial. □ 

Lemma 3. Let[B^,T, P>] be a pre-regular semi-algebraic 
system o/Q[u, y]. One can decide whether its zero set is 
empty or not. If it is not empty, then one can compute a 
regular semi-algebraic system [Q,T, P>] whose zero set in 
R" is the same as that of[B^,T, P>]. 

Proof. If T = 0, we can always test whether the zero 
set of [B^,P>] is empty or not, for instance using CAD. If 
it is empty, we are done. Otherwise, defining Q = B^ A P>, 
the triple [Q, T, P>] is a regular semi-algebraic system. If 
T is not empty, we solve the quantifier elimination problem 
3y(P(u) / 0,T(u,y) = 0,P(u,y) > 0) and let Q be the 
resulting formula. If Q is false, we are done. Otherwise, by 
Lemma[2] above each connected component of B(u) 7^ 0, the 
number of real zeros of the system [B-t,T, P>] is constant. 
Then, the zero set defined by Q is the union of the connected 
components of B(u) 7^ above which [B^,T, P>] possesses 
at least one solution. Thus, Q defines a nonempty open set 
of R d and [Q, T, P>] is a regular semi-algebraic system. □ 

Theorem 1. Let 6 be a semi-algebraic system o/Q[x]. 
Then one can compute a (full) triangular decomposition of 
&, that is, as defined in the introduction, finitely many reg- 
ular semi-algebraic systems such that the union of their zero 
sets is the zero set of 6 . 

Proof. It follows from Lemma [1] and [3] □ 

4. COMPLEXITY RESULTS 

We start this section by stating complexity estimates for 
basic operations on multivariate polynomials. 

Complexity of basic polynomial operations. Let p, q £ 

Q[x] be polynomials with respective total degrees S p ,S q , and 
let x £ x.. Let h p ,h q ,h pq and h r be the height (that is, the 
bit size of the maximum absolute value of the numerator or 
denominator of a coefficient) of p, q, the product pq and the 
resultant res(p, q, x), respectively. In 114] . it is proved that 
gcd(p, q) can be computed within 0(n 2S+1 h 3 ) bit operations 
where S = max(8 p ,S q ) and h = max(fi p , h q ). It is easy to 
establish that h pq and h r are respectively upper bounded by 
h p + h q + nlog(min(5 p , 5 q ) + 1) and 5 q h p + S p h q +nS q \og(8 p + 
1) + nS p log(<5 ? + 1) + log ((S p + S q )\). Finally, let M be a 
k x k matrix over Q[x]. Let S (resp. h) be the maximum 
total degree (resp. height) of a polynomial coefficient of M. 
Then det(M) can be computed within 0(k 2n+5 (S + l) 2n h 2 ) 
bit operations, see [15] . 

We turn now to the main subject of this section, that is, 
complexity estimates for a lazy triangular decomposition of a 
polynomial system under some genericity assumptions. Let 
F C Q[x]. A lazy triangular decomposition (as defined in the 
Introduction) of the semi-algebraic system & = [F, 0,0,0], 



Algorithm 1: LazyRealTriangularize(S) 

Input: a semi-algebraic system & — [F, 0, 0, 0] 
Output: a lazy triangular decomposition of & 

1 X := Triangularize(P) 

2 for Ti el do 

3 bpi := BorderPolynomial(T;, 0) 

4 solve 3y(bpi(u) / 0, Ti(u, y) = 0), and let Qi be the 
resulting quantifier-free formula 
if Q» + false then output [Qi,Ti, 0] 



which only involves equations, is obtained by the above al- 
gorithm. 

Proof of Algorithm [TJ The termination of the algorithm is 
obvious. Let us prove its correctness. Let Pi = [Qi,Ti,%], 
for i = 1 • • • t be the output of Algorithm [T] and let Tj for 
j = t + 1 • • • s be the regular chains such that Qj = false. 
By Lemma [3] each Ri is a regular semi-algebraic system. 
For i = 1 • • • s, define P = sat(Ti). Then we have V(F) = 
L>1 =1 V(Fi), where each Fi is equidimensional. For each i = 
1 • • • s, by Proposition [1] we have 

V(Fi) \ V(b Pi ) = V(Ti) \ V{bpi). 

Moreover, we have 



V{Fi) = (V(Fi) \ V(bpi)) U V(Fi U {bp,}). 



Hence, 



Z R (Ri) = Z R (Ti) \ Z m (b Pi ) C Z R (Fi) C Z m (F) 

holds. In addition, since bpi is regular modulo Fi, we have 

Z m (F)\ui =1 Z R (Ri) = U 1 s = i^m(P) \uLi Z R (Ri) 

C U s i=1 Zu(Fi)\(Zu(Ti)\Zu{bpi)) 
C Ul =l Z M (FiU{bpi}), 

and dim(Uf =1 V(P U {bp,})) < dim(V(F)). So the P i; , for 
i = 1 • • • t, form a lazy triangular decomposition of &. □ 

In this section, under some genericity assumptions for P, 
we establish running time estimates for Algorithm [1] see 
Proposition [3] This is achieved through: 

(1) Proposition [2] giving running time and output size es- 
timates for a Kalkbrener triangular decomposition of 
an algebraic set, and 

(2) Theorem [5] giving running time and output size esti- 
mates for a border polynomial computation. 

Our assumptions for these results are the following: 
(Ho) V(F) is equidimensional of dimension d, 
(Hi) xi, . . . ,Xd are algebraically independent modulo each 
associated prime ideal of the ideal generated by P in 

QM, 

(H2) P consists of m := n — d polynomials, fi, . . . , f m . 
Hypotheses (Ho) and (Hi) are equivalent to the existence of 
regular chains Ti, . . . , T e of Q[a;i, . . . , x n ] such that xi, . . . , Xd 
are free w.r.t. each of Ti , . . . , T e and such that we have 
V(F) = W{T X ) U ■ ■ ■ U W(T e ). 

Denote by S, h respectively the maximum total degree and 
height of /1, .. . ,f m . In her PhD Thesis [22], A. Szanto de- 
scribes an algorithm which computes a Kalkbrener triangu- 
lar decomposition, Ti, . . . , T e , of V(F). Under Hypotheses 

-,0(l)^0(n 2 )-jd+l 



(Ho) to (H2), this algorithm runs in time 



counting operations in Q, while the total degrees of the poly- 
nomials in the output are bounded by n5 0( - m \ In addition, 



Tx,...,T e are square free, strongly normalized [IB] and re- 
duced [T|. 

From Ti , . . . , T e , we obtain regular chains E\, . . . , E B form- 
ing another Kalkbrener triangular decomposition of V(F), 
as follows. Let i = 1 • ■ • e and j = (d + 1) • • • n. Let tij 
be the polynomial of T, with Xj as main variable. Let ei,j 
be the primitive part of ti t j regarded as a polynomial in 
(Q[xi , . . . , X d ] [x d j^\ , • • • , X n ]. Define Ei = {e itd +i, ■ ■ ■ , e i>n }. 
According to the complexity results for polynomial opera- 
tions stated at the beginning of this section, this transfor- 
mation can be done within <5°' m '°' n ' operations in Q. 

Dividing by its initial we obtain a monic polynomial 
di,j of Q(an, . . . ,x d )[x d +i, ...,x„]. Denote by A the regular 
chain {d i} d+i, ■ ■ ■ , rfi.n}- Observe that Di is the reduced lex- 
icographic Grobner basis of the radical ideal it generates in 
Q(xi, . . . , Xd)[xd+i, ■ ■ ■ , x n ]. So Theorem 1 in [12] applies to 
each regular chain Di. For each polynomial di,j, this theo- 
rem provides height and total degree estimates expressed as 
functions of the degree [7] and the height [191 116] of the alge- 
braic set W(Di). Note that the degree and height of W(Di) 
are upper bounded by those of V(F). Write di,j — jj, 

where each /i G <Q[xd+i, ■ ■ ■ , x n ] is a monomial and a M , /3 M 
are in Q[a;i, . . . , x d ] such that gcd(a M , fi^) = 1 holds. Let 7 
be the 1cm of the P^s. Then for 7 and each a^: 

• the total degree is bounded by 25 2m and, 

• the height by 0(S 2m {mh + dm\og(6) + ralog(n))). 
Multiplying di,j by 7 brings e»,j back. We deduce the height 
and total degree estimates for each eij below. 

Proposition 2. The Kalkbrener triangular decomposition 
Ei,...,E e ofV(F) can be computed in «5°(™ 4 )o(n) opem _ 
tions in Q. In addition, every polynomial e%,j has total degree 
upper bounded by 48 2m + 8 m , and has height upper bounded 
by 0{5 2m {mh + dmlog(5) + nlog(n))). 

Next we estimate the running time and output size for 
computing the border polynomial of a regular system. 

Theorem 2. Let R = [T, P] be a squarefree regular sys- 
tem o/Q[u, y], with m — #T and £ = #P. Let bp be the bor- 
der polynomial of R. Denote by S R , h R respectively the max- 
imum total degree and height of a polynomial in R. Then the 
total degree of bp is upper bounded by (£ + m)2 m ^ 1 5 R m , and 
bp can be computed within (n£ + nm) 0{ - n \28 B ) 0{n)0{rn) h R 3 
bit operations. 

Proof. Define G := P U {der(£) | t e T}. We need to 
compute the £+m iterated resultants res(g, T), for all g G G. 
Let g £ G. Observe that the total degree and height of g are 
bounded by 5 R and h R +log(8 R ) respectively. Define r m +i := 
g, n := res(ti,r i+1 ,yi), . . . , n := res(ti, r 2 , J/i). Let i € 
{1, . . . , m}. Denote by Si and hi the total degree and height 
of respectively. Using the complexity estimates stated at 
the beginning of this section, we have 8, < 2 m ~' +1 8 R m ~ l+2 
and hi < 28i+i(hi+i +nlog(<5;+i + 1)). Therefore, we have 

hi < (25fl)° (m 'n° (m) fti{. From these size estimates, one 
can deduce that each resultant r« (thus the iterated resul- 
tants) can be computed within (2fe) 0(mn)+0(m2) n 0(m) fifl 2 
bit operations, by the complexity of computing a determi- 
nant stated at the beginning of this section. 

Hence, the product of all iterated resultants has total 
degree and height bounded by (£ + m)2 m - 1 S R m and (£ + 
m)(25 R ) c "> m2) n 0(m) h R , respectively. Thus, the primitive 



and squarefree part of this product can be computed within 
{nl + nm) 0{n \25 R ) 0{ - n)0{m) h R :i bit operations, based on the 
complexity of a polynomial gcd computation stated at the 
beginning of this section. □ 

Proposition 3. From the Kalkbrener triangular decom- 
position Ei , . . . , E e of Proposition [3 a lazy triangular de- 
composition of fx = • • ■ = f m = can be computed in 

h ^ 1 ' bit operations. Thus, a lazy triangu- 
lar decomposition of this system is computed from the input 
polynomials in singly exponential time w.r.t. n, counting op- 
erations in Q. 

Proof. For each i € {1 • • ■ e}, let bpi be the border poly- 
nomial of [Ei, 0] and let h Ri (resp. 8 Ri ) be the height (resp. 
the total degree) bound of the polynomials in the pre-regular 
semi-algebraic system Ri — [{bpi}^, Ei, 0]. According to Al- 
gorithm [1] the remaining task is to solve the QE problem 
3y(6pi(u) 7^ 0, Ei(u, y) = 0) for each i G {1 • • • e}, which 
can be solved within ((m + l)8 Ri ) 0< ~ dm ^ fi^' 1 ' bit operations, 
based on the results of [20]. The conclusion follows from the 
size estimates in Proposition [2] and Theorem [2] □ 

5. QUANTIFIER ELIMINATION BY RRC 

In the last two sections, we saw that in order to compute 
a triangular decomposition of a semi-algebraic system, a key 
step is to solve the following quantifier elimination problem: 

3y(B(u) + 0,T(u,y) = 0,P(u,y) > 0), (1) 

where [B-t,T, P>] is a pre-regular semi-algebraic system of 
Q[u, y]. This problem is an instance of the so-called real 
root classification (RRC) [27] . In this section, we show how 
to solve this problem when B is what we call a fingerprint 
polynomial set. 

Fingerprint polynomial set. Let R := [B^, T, P>] be a pre- 
regular semi-algebraic system of Q[u, y]. Let D C Q[u]. Let 
dp be the product of all polynomials in D. We call D a 
fingerprint polynomial set (FPS) of R if: 
(i) for all a G K d , for all b G B we have: 

dp{a) ± =^> b(a) 0, 
(ti) for all a, /3 G R d with a / /3 and dp(a) / 0, dp(/3) / 0, 
if the signs of p(a) and p(/3) are the same for all p G D, 
then R(a) has real solutions if and only if R(/3) does. 
Hereafter, we present a method to construct an FPS based 
on projection operators of CAD. 

Open projection operator 1211 [4]. Hereafter in this sec- 
tion, let u = tti < • ■ • < Ud be ordered variables. Let 
p G Q[u] be non-constant. Denote by factor(p) the set of 
the non-constant irreducible factors of p. For A C Q[u], 
define factor(A) = U p gA factor(p). Let C d (resp. Co) be 
the set of the polynomials in factor(p) with main variable 
equal to (resp. less than) u d . The open projection operator 
(oproj) w.r.t. variable Ud maps p to a set of polynomials of 
Q[ui, . . . , Ud-\] defined below: 

oproj(p,ti d ) := Co U U/, 9G c d , f* g factor(res(/, g, u d )) 
U U/gc d factor(init(/, Ud) ■ discrim(/, Ud))- 

Then, we define oproj(A, u d ) := oproj (n pg yt p, u d ). 

Augmentation. Let A C Q[u] and x G {ui, ■ ■ ■ , u d }. Denote 
by der(^4, x) the derivative closure of A w.r.t. x, that is, 
der(^,a;) := U peA {der w (p, x) | < i < deg(p, x)}. The 



open augmented projected factors of A is denoted by oaf (A) 
and denned as follows. Let k be the smallest positive integer 
such that A C Q[ui, . . . , Ufc] holds. Denote by C the set 
factor(der(A, Ufc)); we have 

• if fc = 1, then oaf (A) := C; 

• if k > 1, then oaf (A) := C U oaf (oproj(C, Ufe)). 

Theorem 3. £et A C Q[u] &e /mite and let a be a map 
from oaf (A) to the set of signs { — 1, +1}- Then the set Sd := 
n p( =oaf(A) {it G R d | p{u) cr(p) > 0} is either empty or a 
connected open set in R d . 

Proof. By induction on d. When d = 1, the conclusion 
follows from Thorn's Lemma [2]. Assume d > 1. If d is not 
the smallest positive integer fc such that A C Q[ui, • • • , Ufc] 
holds, then Sd can be written Sd-i x R and the conclusion fol- 
lows by induction. Otherwise, write oaf (A) as C U E, where 
C — factor(der(A, Ud)) and E = oaf(oproj(C, Udj)- We have: 
-E C Q[«i,-- - ,Ud-i]- Denote by M the set (l p eE {u G 
R''" 1 | p(u)a(p) > 0}. If M is empty then so is S d and 
the conclusion is clear. From now on assume M not empty. 
Then, by induction hypothesis, M is a connected open set 
in R d_1 . By the definition of the operator oproj, the prod- 
uct of the polynomials in C is delineable over M w.r.t. Ud- 
Moreover, C is derivative closed (may be empty) w.r.t. Ud- 
Therefore n peoaf( A) {u G R d | p(u) a(p) > 0} C M X R is ei- 
ther empty or a connected open set by Thorn's Lemma. □ 

Theorem 4. Let R := [B^,T, P>] be a pre-regular semi- 
algebraic system o/Q[u, y]. The polynomial set oaS(B) is a 
fingerprint polynomial set of R. 

PROOF. Recall that the border polynomial bp of [T, P] di- 
vides the product of the polynomials in B. We have factor(_B) 
C oaf(B). So oaf(B) satisfies (i) in the definition of FPS. 
Let us prove (ii). Let dp be the product of the polynomi- 
als in oaf(B). Let a, /3 G R d such that both dp(a) / 0, 
dp(/3) hold and the signs of p(a) and p(/3) are equal for 
all p G oaf(B). Then, by Theorem [3l a and /3 belong to 
the same connected component of dp(u) 0, and thus to 
the same connected component of B(u) 0. Therefore the 
number of real solutions of R(a) and that of R(fi) are the 
same by Lemma [5] □ 

From now on, let us assume that the set B in the pre- 
regular semi-algebraic system R — [Bjt,T, P>] is an FPS 
of R. We solve the quantifier elimination problem (fTJ) in 
three steps: (si) compute at least one sample point in each 
connected component of the semi-algebraic set defined by 
B(u) 7^ 0; (S2) for each sample point a such that the spe- 
cialized system R(ct) possesses real solutions, compute the 
sign of 6(a) for each b G B; (S3) generate the corresponding 
quantifier-free formulas. 

In practice, when the set B is not an FPS, one adds some 
polynomials from oaf(B), using a heuristic procedure (for 
instance one by one) until Property (ii) of the definition of an 
FPS is satisfied. This strategy is implemented in Algorithm^ 
of Section [6] 

6. IMPLEMENTATION 

In this section, we present algorithms for LazyRealTrian- 
gularize and RealTriangularize that we have implemented on 
top of the Regular/Chains library in Maple. We also provide 
experimental results for test problems which are available at 
www. orcca. on. ca/~cchen/issaclO . txt 



Algorithm 2: GeneratePreRegularSas(S) 

Input: a semi-algebraic system & — [F, N>, P > , H^] 
Output: a set of pre-regular semi-algebraic systems 



[Bi ^i,Ti, Pi~ 
Z R (<3) = 



i = l. 



. e, such that 



Uf = l^K(-Bi^, Ti, Piy) 

Uf =1 Z K (sat(T i ) U {U beBz b},N>,P > ,H 7 ,) 

1 X := Triangularize(F); X' := 

2 for p G P U H do 



for T G X do 

for C G Regularize^, T) do 
|_ if p £ sat(C) then X' := X' U {C} 

X:=X'; X' := 



7 X:= {[T,0] I T G X}; X' := 

8 for p G N do 



9 
10 
11 
12 

13 
14 



for [T, N'} G X do 

for C G Regularize^, T) do 
if p G sat(C) then 
I %' :=%'U{[C,N'}} 
gIsg 

L X':=X'U{[C,7V'U{ P }]} 



X:=X'; X' := 

16 X := {[T, N', P, H] I [T, N'] G X} 

17 for [T,N',P,H] G X do 



18 
19 



B := BorderPolynomialSet(T, N' U P U H) 
output [B, T, N' U P] 



Algorithm 3: GenerateRegularSas(B, T, P) 

Input: 6 = [B^,T, P>], a pre-regular semi-algebraic 
system of Q[u, y], where u = ui, . . . , Ud and 

y = yi, ■■ -,yn-d- 

Output: A pair (D,1Z) satisfying: 

(1) D C Q[u] such that factor(B) C D; 

(2) 1Z is a finite set of regular semi-algebraic systems, 
s.t. yj Relz Zw.{R) = ^ K ( J D # ,T,P>). 

1 D := factor(B \ Q) 

2 if d = then 
if RealRootCounting(T,P) = then 

I return (D, 0) 
else 

[_ return (D, {[true, T, P]}) 



7 while true do 

8 S:= SamplePoints(D,d); Go := 0; Gi := 

9 for s G S do 

10 if RealRootCounting(T(s),P(s)) = then 

11 J Go := Go U {GenerateFormula(Z), s)} 

12 else 

13 |_ Gi := Gi U {GenerateFormula(Z), s)} 

14 if Go n Gi = then 

15 Q := Disjunction(Gi) 

16 if Q = /a/se then return (D, 0) 

17 else return (D, {[Q, T, P]}) 

18 else 

19 select a subset D' C oaf (B) \ D by some 
heuristic method 
D := D U D' 



Basic subroutines. For a zero-dimensional squarefree regu- 
lar system [T,P], RealRootCounting(T, P) [23] returns the 
number of real zeros of [T, P>]. For A C Q[t*i, • • • , uj] 
and a point s of Q d such that p(s) 7^ for all p G A, 
GenerateFormula(A, s) computes a formula A p gA (pff P , s >0), 
where ct PiS is defined as+lifp(s)>0 and —1 otherwise. For 
a set of formulas G, Disjunction(G) computes a logic formula 
$ equivalent to the disjunction of the formulas in G. 

Proof of Algorithm [2] Its termination is obvious. Let us 
prove its correctness. By the specification of Triangularize 
and Regularize, at line 16, we have 

Z{F,P^ U H±) = U IT ,Ar',P,H]6i^(Bat(T), P^ U 

Write V-)[t,n' ,p,h]<=% as Ut- Then we deduce that 

Z R (F,N>,P > ,H^) = U T Z u (sat(T),N>,P > ,H 9t ). 

For each [T, N' , P, H] , at line 19, we generate a pre-regular 
semi-algebraic system [B^,T, iV> UP>]. By Proposition [1] 
we have 

Z M (sat(T),iV>,P > , J H'^) = 

Z R (B^, T, N'> U P>) U Z R (sat(T) U {H beB b}, N>,P > ,H^), 
which implies that 



Algorithm 6: RealTriangularize(S) 



So Algorithm [2] satisfies its specification 



U T Z R (B^,T,N' > UP > ) 
U T Z R (sat(r) U {IWb}, A/>, P>, P^) 



Algorithm 4: SamplePoints(A, k) 



Input: A C Q[jci, • ■ • , Xk] is a finite set of non-zero 

polynomials 
Output: A finite subset of Q fe contained in 
(IIpgAp) 7^ and having a non-empty intersection with 
each connected component of (H p ^ap) 7^ 0. 
1 if k = 1 then 

return one rational point from each connected 
component of H P £A p 7^ 

3 else 

4 Ak := {p e A I mvar(p) = x fe }; A' := oproj(yl, x k ) 

5 for s G SamplePoints(A', k — 1) do 

6 Collect in a set S one rational point from each 
connected component of H pe A k p{s, x k ) 7^ 0; 
for ft £ S do output (s, a) 



Algorithm 5: LazyRealTriangularize(S) 

Input: a semi-algebraic system 6 = [F, N> , P>, H-t] 
Output: a lazy triangular decomposition of © 

1 T := GeneratePreRegularSas(F,iV,P,P) 

2 for [B,T,P'] G I do 

3 (D, Tl) = GenerateRegularSas(B, T, P') 

4 if 72. 7^ then output TZ 



Proof of Algorithms [3] and 131 By the definition of oproj, 
Algorithm [4] terminates and satisfies its specification. By 
Theorem [H oaf(S) is an FPS. Thus, by the definition of an 
FPS, Algorithm [3] terminates and satisfies its specification. 

Proof of Algorithm [5] Its termination is obvious. Let us 
prove the algorithm is correct. Let Ri, i = 1 ■ ■ ■ t be the 



Input: a semi-algebraic system (5 = [F, N> , P> , H^] 
Output: a triangular decomposition of S 

1 X := GeneratePreRegularSas(F,iV,P,P) 

2 for [B, T, P'\ G X do 

3 (D, TZ) = GenerateRegularSas(B, T, P') 

4 if 72. 7^ then output 

5 for p £ D do 

6 [_ output RealTriangularize(P U {p}, AT, P, 77) 



output. By the specification of each sub-algorithm, each Ri 
is a regular semi-algebraic system and we have: 

ut =1 Z u (Ri) C Z„(6). 

Next we show that there exists an ideal I C Q[x], whose 
dimension is less than dim(.Z(F, P^ U H^)) and such that 
Z M (6) \ U* =1 Z R (P0 C Z m (T) holds. 

At line 1, by the specification of Algorithm [2] we have 

Z R (6) = U T Z R (B # ,T,P^) 

U T ZB(sat(T) U {n 6es &}, iV>,P>,P^). 

At line 3, by the specification of Algorithm^ for each B, we 
compute a set D such that factor(P) C D and 

U T Z B (Py, T,P>) = U* = A(Pi) 

both hold. Combining the two relations together, we have 

Z R (S) = U T Zu(Ri) 

U T ^n(sat(r) U {Yl peD p}, N> , P> , P # ). 

Therefore, the following relations hold 

Z R (6)\uUiZR(R t ) 
C U T ^M(sat(T)u{n p6D p},Af>,P > ,P^) 
C Z R (n T (sat(T) u {n peD p})) . 

Define 

j = n T (sat(r)u{n pSD p}). 

Since each p G B is regular modulo sat(T), we have 
dim(I) < dim (n T sat(T)) < dim(Z(P, P^ U B^)). 

So all Ri form a lazy triangular decomposition of &. □ 

Proof of Algorithm [6] For its termination, it is sufficient 
to prove that there are only finitely many recursive calls 
to RealTriangularize. Indeed, if [F, N, P, H] is the input of a 
call to RealTriangularize then each of the immediate recursive 
calls takes [f U{p}, N, P, H] as input, where p belongs to the 
set D of some pre-regular semi-algebraic system [D^, T, P>]. 
Since p is regular (and non-zero) modulo sat(T) we have: 

<F> C <FUM>. 

Therefore, the algorithm terminates by the ascending chain 
condition on ideals of Q[x]. The correctness of Algorithm [6] 
follows from the specifications of the sub-algorithms. □ 

Table 1. Table 1 summarizes the notations used in Tables 2 
and 3. Tables 2 and 3 demonstrate benchmarks running in 
Maple 14 /3 1, using an Intel Core 2 Quad CPU (2.40GHz) 
with 3.0GB memory. The timings are in seconds and the 
time-out is 1 hour. 

Table 2. The systems in this group involve equations only. 
We report the running times for a triangular decomposition 



Table 1 Notations for Tables 2 and 3 



symbol 


meaning 




mi m r»r>r f\T r»nn ^tiniiQ in tnf 1 mtinf evctpm 

llLllllULi Ul 1_.U LldHjlUllO 111 Ull\_. IIIULIL 3 V QLljlll 


7T v 


number of variables in the input equations 


d 


maximum total degree of an input equation 


G 


Groebner:-Basis (plex order) in Maple 


T 


Triangularize in RegularChains library of Maple 


LR 


LazyRealTriangularize implemented in Maple 


R 


RealTriangularize implemented in Maple 


Q 


Qepcad b 


> lh 


computation does not complete within 1 hour 


FAIL 


Qepcad b failed due to prime list exhausted 



Table 2 Timings for varieties 



system 


#v/#e/d 


G 


T 


LR 


Hairer-2-BGK 


13/ 11/ 4 


25 


1.924 


2.396 


Collins-jsc02 


5/ 4/ 3 


876 


0.296 


0.820 


Leykin-1 


8/ 6/ 4 


103 


3.684 


3.924 


8-3-config-Li 


12/ 7/ 2 


109 


5.440 


6.360 


Lichtblau 


3/ 2/ 11 


126 


1.548 


11 


Cinquin-3-3 


4/ 3/ 4 


64 


0.744 


2.016 


Cinquin-3-4 


4/ 3/ 5 


> lh 


10 


22 


DonatiTraverso-rev 


4/ 3/ 8 


154 


7.100 


7.548 


Cheaters-homotopy- 1 


7/ 3/ 7 


3527 


174 


> lh 


hereman-8.8 


8/ 6/ 6 


> lh 


33 


62 


L 


12/ 4/ 3 


> lh 


0.468 


0.676 


dgp6 


17/19/ 2 


27 


60 


63 


dgp29 


5/ 4/ 15 


84 


0.008 


0.016 



of the input algebraic variety and a lazy triangular decom- 
position of the corresponding real variety. These illustrate 
the good performance of our tool. 

Table 3. The examples in this table are quantifier elimi- 
nation problems and most of them involve both equations 
and inequalities. We provide the timings for computing a 
lazy and a full triangular decomposition of the correspond- 
ing semi-algebraic system and the timings for solving the 
quantifier elimination problem via Qepcad b [5] (in non- 
interactive mode). Computations complete with our tool on 
more examples than with Qepcad b. 

Remark. The output of our tools is a set of regular semi- 
algebraic systems, which is different than that of Qepcad B. 
We note also that our tool is more effective for systems with 
more equations than inequalities. 

Acknowledgments. The authors would like to thank the 
referees for their valuable remarks that helped to improve 
the presentation of the work. 



Table 3 Timings for semi-algebraic systems 



system 


#v/#c 


/d 


T 


LR 


R 


Q 


BM05-1 


4/ 2/ 


3 


0.008 


0.208 


0.568 


86 


BM05-2 


4/ 2/ 


4 


0.040 


2.284 


> lh 


FAIL 


Solotareff-4b 


5/ 4/ 


3 


0.640 


2.248 


924 


> lh 


Solotareff-4a 


5/ 4/ 


3 


0.424 


1.228 


8.216 


FAIL 


putnam 


6/ 4/ 


2 


0.044 


0.108 


0.948 


> lh 


MPV89 


6/ 3/ 


4 


0.016 


0.496 


2.544 


> lh 


IBVP 


8/ 5/ 


2 


0.272 


0.560 


12 


> lh 


Lafferriere37 


3/ 3/ 


4 


0.056 


0.184 


0.180 


10 


Xia 


6/ 3/ 


4 


0.164 


191 


739 


> lh 


SEIT 


11/ 4/3 


0.400 


> lh 


> lh 


> lh 


p3p-isosceles 


7/ 3/ 


3 


1.348 


> lh 


> lh 


> lh 


p3p 


8/ 3/ 


3 


210 


> lh 


> lh 


FAIL 


Ellipse 


6/ 1/ 


3 


0.012 


> lh 


> lh 


> lh 
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